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We report on a CDF measurement of the total cross section and rapidity distribution, da/dy, 
for qq — > 7*/^ — > e"'"e~ events in the Z boson mass region (66 < Mee < 116 GeV/c^) produced 
in pp collisions at ^/s = 1.96 TcV with 2.1 fb~^ of integrated luminosity. The measured cross 
section of 257± 16 pb and da/dy distribution are compared with Next-to-Leading-Order(NLO) and 
Ncxt-to-Ncxt-to-Lcading-Ordcr(NNLO) QCD theory predictions with CTEQ and MRST/MSTW 
parton distribution functions (PDFs). There is good agreement between the experimental total 
cross section and da/dy measurements with theoretical calculations with the most recent NNLO 
PDFs. 
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I. INTRODUCTION 

Accurate predictions using perturbative quantum 
chromodynamics (QCD) are critical for understanding 
experimental results at hadron colliders. Such predic- 
tions depend on the accuracy of input parton distribu- 
tion functions (PDFs), which at present cannot be cal- 
culated and are obtained from analysis of data from a 
broad range of processes. Precise knowledge of PDFs 
will be particularly important for analysis of data at 
the Large Hadron Collider (LHC) where new phenom- 
ena may be revealed via small deviations from Standard 
Model (SM) predictions. The Drell-Yan process in 
which quark-antiquark annihilations form intermediate 
7* or Z {"f* /Z) vector bosons decaying to lepton pairs, is 
particularly useful in providing information on PDFs at 
— M|^, where Mu is the invariant mass of the dilep- 
ton pair. In the leading order (LO) approximation, the 
momentum fractions xi, X2 carried by the initial state 
quarks in the proton and antiproton, respectively, are 
related to the rapidity y [1] of the 7*/^ boson via the 
equation xi^2 = {Mee/ ^/s)e^''' , where -^s is the center of 
mass energy. Dilepton pairs produced at large y origi- 
nate from collisions in which one parton carries a large 
and the other a small momentum fraction x. A mea- 
surement of da/dy at large y tests PDFs at high x, a 
region not well constrained by current results. There- 
fore, precise measurements of W and Z boson rapidity 
distributions at the Tevatron determine the size of higher 
order QCD terms and can be used to further refine cur- 
rent PDF models. Furthermore since the Z production 
cross section is predicted with an accuracy of ~ 2% Q, 
precise measurements of the rate of Z production at the 
Tevatron and the LHC can be used to determine the inte- 
grated luminosity Q more precisely than the traditional 
method of using the total inelastic cross section. This 
has particular applicability to sub-processes initiated by 
a quark and an anti-quark and can reduce the uncertainty 
in the determination of LHC and Tevatron cross sections. 

The most recent Tevatron measurement of da/dy for 
e~^e^ pairs in the Z boson mass region was performed 
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by the DO Q experiment, using a data-set correspond- 
ing to 0.4 fb"'^ of integrated luminosity. Here, we re- 
port on a new measurement of da/dy at the Tevatron 
with an integrated luminosity of 2.1 fb~^. The mea- 
sured rapidity range extends to \y\ ~ 2.9, close to the 
kinematic limit of \y\ — 3.0 for Z boson production at 
^/s = 1.96 TeV. The da/dy distribution is compared 
to the predictions of perturbative QCD calculations 
in Next-to-Leading-Order (NLO) and Next-to-Next-to- 
Leading-Order (NNLO) with different PDF models. 



II. EVENT SELECTION AND ANALYSIS 
METHOD 

The data sample corresponds to an integrated lu- 
minosity of 2.1 fb^^ collected by the CDF H Detec- 
tor at Fermilab [Sl during 2004-2007. CDF H uses a 
1.4 T solenoidal magnetic spectrometer surrounded by 
projective-tower-geometry calorimeters and outer muon 
detectors. Charged particle directions and momenta are 
measured by an open-cell drift chamber (COT), a sili- 
con vertex detector (SVX), and an intermediate silicon 
layer (ISL). The coverage of COT tracking in pseudo- 
rapidity is |r?| < 1.2 Q. Reconstructed tracks are used 
to determine the pp collision point along the beam line 
(z„ertex), which is required to be within z = ±60 cm of 
the detector. The energies and directions [2] of electrons, 
photons, and jets are measured by two separate calorime- 
ters: central (|?7| < 1.1) and plug (1.1 < |?7| < 3.6). Each 
calorimeter has an electromagnetic (EM) compartment 
with a shower maximum detector followed by a hadronic 
(HAD) compartment. Three topologies of e'^e" pairs 
are considered: two central electrons (CC), one central 
and one plug electron (CP), and two plug electrons (PP). 
The inclusion of PP events allows the measurement of Z 
bosons in the forward rapidity region which corresponds 
to high and low parton momentum fractions. 

Data are collected using a three-level trigger system 
,6;] and trigger paths with either one central electron 
or two electrons (central or plug) with transverse en- 
ergy Et > 18 GeV. Electron identification requirements 
are imposed to select signal events and to suppress 
background. Both electron candidates are required to 
be isolated from any other calorimetric activity. The 
fraction of energy in the HAD calorimeter towers behind 
the EM shower is required to be small 01 j as expected 
for an EM shower. Electron candidates with Et > 25 
GeV for CC and PP events, and Et > 20 GeV for CP 
events, are selected in the central (|7/| < 1.1), and plug 
(1.2 < \ri\ < 2.8) fiducial regions of the calorimeters. 
Central electron candidates must have a COT track that 
extrapolates to a shower cluster in the EM calorimeter 
and a track momentum consistent with the calorimeter 
measurement. Central and plug electron candidates are 
required to have EM-like transverse shower profiles us- 
ing the shower maximum detectors. In order to reduce 
background we require that at least one of the plug elec- 
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trons in PP events has a track reconstructed in the SVX 
that points to the EM cluster in the calorimeter. The effi- 
ciency of having at least one electron matched to an SVX 
track is about 85%. The selected number of CC, CP, and 
PP events with 66 < M^e < 116 GeV/c^ is 50752, 86203, 
and 31415, respectively. 



using the PYTHIA [8| generator combined with a GEANT 



ll| simulation of the CDF detector. 
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FIG. 1: The distribution of background events as a function 
of rapidity for CC, CP and PP e^e~ candidates (shown as a 
fraction of all candidate events selected in this analysis). The 
error bars include the statistical and systematic uncertainty. 



III. BACKGROUNDS 

The main backgrounds are QCD dijet and photon plus 
jet events in the plug region (because of the limited track- 
ing at large \r]\). The jet background is measured sepa- 
rately in each e+e^ pair topology by statistically sepa- 
rating electrons from jets on the basis of the transverse 
energy profile distributions in the calorimeter and the 
invariant mass distribution of e+e^ pairs Q- The distri- 
bution of background events as a function of rapidity for 
CC, CP and PP e^e~ candidates (shown as a fraction 
of all candidate events selected in this analysis) is shown 
in Fig. [H The fractional contribution of the total back- 
ground to the number of selected events is 0.24 ± 0.03% 
(stat ® syst) for CC, 1.55±0.44% for CP, and 3.40±0.75% 
for PP events. The background from electroweak {WW, 
WZ, 14^-l-jets, and Z — > r+r^) and tt processes is esti- 
mated from simulation to be 0.41 ± 0.02%. 



IV. ACCEPTANCE AND EFFICIENCIES 

The acceptance is defined as the ratio of the number of 
Monte Carlo (MC) simulation events that pass selection 
criteria in each y bin of the reconstructed final state e^e~ 
pair (including resolution smearing) to the number of MC 
generated events in each true y bin of the generated /Z 
boson. The resolution in the measurement of the e"*"e^ 
invariant mass is 2.2 GeV/c^, and the resolution in the 
measurement of y is 0.015. The acceptance is modeled 



The PYTHIA generator includes a LO QCD interac- 
tion (g + 5 — > "f*/Z), initial state QCD radiation, parton 
shower fragmentation, the j* /Z — >■ e'^e" decay, and pho- 
ton radiation from the final state. The version of pythia 
used at CDF has additional ad-hoc tuning in order to 
accurately represent the Z boson transverse momentum 
distribution measured in data. To reconstruct the sim- 
ulated events in the same way as data, the calorimetry 
energy scale, resolutions, and selection efficiencies used 
in the detector simulation are tuned using data. 

Because the acceptance depends on modeling of the Z 
boson rapidity, transverse momentum and angular dis- 
tributions of the electron pairs, it is important to correct 
for possible model dependences arising from the choice 
of the event generator or a particular PDF set. The un- 
corrected acceptance is calculated using the CTEQ5L 3 
LO PDFs, and we compare relevant kinematic distribu- 
tions in the MC simulation to those observed in the data 
to correct the acceptance for possible observed discrep- 
ancies. 

While the generated rapidity spectrum is in good 
agreement with the data for y < 2, the data and sim- 
ulation do not agree at larger values of y. To correct 
for this discrepancy, we modify the MC generated event 
spectrum (dN/dy) so that the final accepted MC spec- 
trum matches the spectrum in data, as shown in Fig. 
[51 A comparison of the reconstructed transverse momen- 
tum spectra of the e'^e~ pairs in the data and the MC 
simulation reveals good agreement as shown in Fig. |31 
Modifying the Pt spectrum in simulation to exactly fol- 
low the data leads to a negligible change in the calculated 
acceptance. A comparison of the average Pt of events in 
bins of y shows that the data is well modeled by the sim- 
ulation. Similarly, a study of the angular distribution 
in 9 (where 9 is the polar angle of the final state elec- 
tron in the CoUins-Soper frame .10]) in data shows good 
agreement with the simulation as illustrated in Fig. 2] 

The acceptance (A) and efficiencies (e) are determined 
as a function of boson rapidity. The contributions of each 
topology to the product Axe arc shown in Fig. [5l 



V. DIFFERENTIAL AND TOTAL CROSS 
SECTIONS 

The differential cross section is given by 



dah*/Z) 
dy 



iy) 



c{y)Aye^ytx^^[{A, x e^(y))ej„g(y)/: 



stgiy) ~ Nbkg{y) is the number of events after 



where N. 

subtracting background, c{y) is a correction factor used 
in order to yield da/dy at the center of the bin, and Ay 
is the y bin size {Ay — 0.1 up to y = 2.7 and Ay = 0.2 
for the last bin, 2.7 < y < 2.9). The sum index i runs 
over the e+e" topologies (CC, CP, PP), Ai x €i{y) is 
the combined acceptance and event selection efficiency. 
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FIG. 2: The ratio of the rapidity distributions dN/ dy of recon- 
structed e^e~ pairs in data to the results from tire simulation 
using PYTHIA, CTEQ5L LO PDFs and CDF W/Z tuning pa- 
rameters @] . (A) PYTHIA before modification of the rapidity 
distribution. (B) pythia after modification of the rapidity 
distribution . 
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FIG. 3: The reconstructed transverse momentum distribution 
of e^e" pairs (after all event selection cuts) compared to the 
results from the simulation using pythia, CTEQ5L LO PDFs 
and CDF W/Z tuning parameters Q- The black points are 
data and the blue solid line is the pythia MC prediction. 




elrigiy) is the trigger efficiency, Ci is the total integrated 
luminosity for each topology, and ezvtx is the acceptance 
for the pp collision vertex to occur within z = ±60 cm of 
the center of the detector. The ezvtx in the data taken 
before June 2006 is 95.8 ± 0.2% and after that is 96.8 ± 
0.2%. 

Systematic uncertainties in da/dy originate from un- 
certainties in the estimates of the acceptance, back- 
grounds, electron identification efficiency, SVX tracking 
efficiency, and modeling of material in the detector. Un- 
certainties associated with correcting the acceptance for 
differences between kinematic distributions in data and 
simulation are found to be negligible. The total system- 
atic uncertainty is ~ 1.0% of da/dy for \y\ < 2.5, in- 
creasing to 10.0% at \y\ = 2.9. The uncertainty on the 
integrated luminosity (lum.) is 6%. 



VI. RESULTS 

The measured da/dy values, which are symmetric 
about y = 0, are shown versus \y\, with statistical and 
systematic uncertainties, in Fig. [B] and Table HI The 



FIG. 4: The reconstructed polar angle distribution of e e~ 
pairs (after all event selection cuts) compared to the results 
from the simulation using pythia, CTEQ5L LO PDFs and 
CDF W/Z tuning parameters Q]. The black points are data 
and the blue solid line is the pythia MC prediction. 



total cross section, derived from integrating da/dy up to 
\y\ = 2.9, is cr = 256.6±0.7(stat.)±2.0(syst.)±15.4(lum.) 
pb. These results are compared to QCD predictions 
at LO with CTEQ5L at NLO ^ with MRST2001 
(NLO) [M], MRST2004 (NLO) [H, CTEQ6.1M (NLO) 
[111, CTEQ6.6M (NLOml, and MSTW2008 (NLO) [| 
PDFs, and at NNLO flsTwith MRST2006E (NNLO) 
[Hi and MSTW2008 (NNLO) 3] PDFs. The measured 
total cross section is consistent with both NLO and 
NNLO calculations as shown in Table [Til 

In comparing the shape of the measured da/dy to the- 
ory, the latter distributions are normalized to the mea- 
sured total cross section of 256.6 pb. The ratios of the 
measured da/dy to the QCD calculations at LO, NLO 
and NNLO with the above mentioned PDFs are shown 
in Figures [71 [51 and [S] The yellow bands in the fig- 
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FIG. 5: The product of kinematic acceptance and event se- 
lection efficiency vs. the rapidity of the e'*'e~ pair. The black 
points are the sum of all topologies. 
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FIG. 6: The measured da/dy for pp — >• Z^'/j* -> e+e~ 
over the entire rapidity range. The points are the mea- 
sured cross section versus |y| and the solid line is the the- 
ory prediction (scaled to the measured total cross section) for 
CTEQ6.1M(NLO) PDFs. 



ures correspond to the uncertainties associated to the 
MSTW2008E (NLO and NNLO) PDFs, which are given 
with 68% C.L. errors and to the those associated to the 
other sets of PDFs, given with 90% C.L. errors. A 
comparison (including statistical and systematic uncer- 
tainties) is shown in Table IIIIl Better agreement is ob- 
tained for the MSTW (2008) PDF at NNLO compared 
to NLO. The NLO CTEQ6.1M, CTEQ6.6M, the NNLO 
MRST2006 and MSTW2008 PDFs all describe the data 
weU. The older NLO MRST (2004) set provides a poorer 
description of the data and, as expected, so does the LO 
PDF, CTEQ5L. The MSTW(2008) PDF used a prelimi- 
nary version of the data, presented in this paper in their 
fit. The correlations 0, [3] between the uncertainties in 
different y bins are included in the comparison. 

In summary, high-statistics measurement of 7*/Z pro- 
duction in the Z mass region and of its rapidity dis- 
tribution are found to agree with theoretical calcula- 
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FIG. 7: The ratio of the experimental distribution of da/dy 
(statistical and systematic uncertainties combined) to the 
theoretical predictions with CTEQ NLO PDFs (CTEQ6.1M 
and CTEQ6.6M). The yellow bands corresponds to the 
CTEQ6.6M PDFs 90% C.L. uncertainties. The test in- 
cludes the data statistical and systematical uncertainties. 



o 

t 

(0 
(0 

a 



• 

o 


NLO MRST2004 {x^ = 70/27, CL=2.05x 10"') 
NLO MSTW2008E (x^ = 51/27, CL=0.005) 
MSTW2008E 68% C.L. PDFs Uncertainties 






































r*-*-rT 


















" , , , , 


, , , 


, , , 




, , , 









0.5 1 1.5 2 2.5 
e^e' pair rapidity 



FIG. 8: The ratio of the experimental distribution of da/dy 
(statistical and systematic uncertainties combined) to the the- 
oretical predictions with MRST2004 and MSTW2008E NLO 
PDFs. The yellow bands corresponds to the MSTW2008E 
PDFs 68% C.L. uncertainties. The test includes the data 
statistical and systematical uncertainties. 



tions that use recent NLO and NNLO PDFs. (A pre- 
liminary version of these results has been used in the 
determination of the most recent PDFs). The pre- 
cise measurement of the total production cross section, 
a = 256.6 ± 0.7(stat.) ± 2.0(syst.) ± 15.4(lum.) pb can be 
used to set the normalization of other processes at the 
Tevatron and the LHC. 
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FIG. 9: The ratio of the experimental distribution of da/dy 
(statistical and systematic uncertainties combined) to the 
theoretical predictions with MRST2006E and MSTW2008E 
NNLO PDFs. The yellow bands corresponds to the 
MSTW2008E PDFs 68% C.L. uncertainties. The test in- 
cludes the data statistical and systematical uncertainties. 



TABLE L Differential cross sections for production of e^e~ 
pairs in the mass range 66 < Mee < 116 GeV/c^. The first 
and second uncertainties are statistical and systematical, re- 
spectively. The 6% luminosity uncertainty is not included. 
The quoted y values correspond to the center of the bin. The 
bin size is 0.1 up to y = 2.7 and 0.2 for the last bin. 



y 


da/dy[pb] 


y 


da / dy [ph] 


0.05 


69.46 ± 0.73 ± 0.49 


1.55 


50.07 ±0.62 ±0.37 


0.15 


71.03 ±0.74 ±0.49 


1.65 


46.59 ±0.61 ±0.35 


0.25 


71.10 ±0.74 ±0.49 


1.75 


40.97 ±0.58 ±0.34 


0.35 


70.01 ±0.72 ±0.48 


1.85 


37.04 ± 0.56 ± 0.33 


0.45 


67.97 ±0.70 ±0.47 


1.95 


33.02 ±0.55 ±0.31 


0.55 


68.22 ± 0.70 ± 0.47 


2.05 


27.65 ±0.52 ±0.25 


0.65 


66.58 ± 0.69 ± 0.47 


2.15 


21.84 ±0.49 ±0.23 


0.75 


66.81 ±0.70 ±0.48 


2.25 


18.35 ±0.50 ±0.20 


0.85 


65.05 ± 0.69 ± 0.49 


2.35 


14.13 ±0.49 ±0.17 


0.95 


64.70 ± 0.69 ± 0.50 


2.45 


8.80 ±0.45 ±0.10 


1.05 


62.74 ±0.67 ±0.50 


2.55 


5.68 ±0.44 ±0.09 


1.15 


62.02 ± 0.66 ± 0.49 


2.65 


2.93 ±0.41 ±0.15 


1.25 


58.80 ± 0.65 ± 0.48 


2.80 


0.87 ±0.22 ±0.11 


1.35 


56.02 ±0.65 ±0.43 


2.95 




1.45 


53.37 ±0.63 ±0.40 
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x'/DOF 


CL 


CTEQ5L(L0) 
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